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We combine results from searches by the CDF and DO Collaborations for a standard model Higgs 
boson (H) in the processes gg — > H — > W + W~ and gg — > H — > ZZ in pp collisions at the Fermilab 
Tevatron Collider at yfs = 1.96 TeV. With 8.2 fb -1 of integrated luminosity analyzed at CDF 
and 8.1 fb" 1 at DO, the 95% C.L. upper limit on a(gg -> H) x B(H -> W+W") is 1.01 pb at 
mjj = 120 GeV, 0.40 pb at ran = 165 GeV, and 0.47 pb at mn = 200 GeV. Assuming the presence 
of a fourth sequential generation of fermions with large masses, we exclude at the 95% Confidence 
Level a standard-model-like Higgs boson with a mass between 124 and 286 GeV. 



Preliminary Results 



* The Tevatron New-Phenomena and Higgs Working Group can be contacted at TEVNPHWG@fnal.gov. More information can be found 
at http://tevnphwg.fnal.gov/ 
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I. INTRODUCTION 

Exploring the mechanism for the breaking of the SU(2) x f (1) electroweak gauge symmetry is a priority in high 
energy physics. Not only are this symmetry and its breaking jl| necessary components for the consistency of the 
successful standard model (SM) [2j, but measurable properties of the breaking mechanism are also very sensitive to 
possible phenomena that have not yet been observed at collider experiments. Measuring these properties, or setting 
limits on them, can constrain broad classes of extensions to the SM. 

A natural extension to the SM that can be tested with Higgs boson search results at the Fermilab Tevatron Collider 
is the addition of a fourth generation of fermions with masses much larger than those of the three known generations [H . 
While fits to precision electroweak data favor a low-mass Higgs boson in the SM, the addition of a fourth generation 
of fermions to the SM modifies the fit parameters such that a much heavier Higgs boson is allowed 0, HI • Currently 
available data are consistent with a Higgs boson of mass 900 GeV with suitably chosen fourth generation fermion 
masses [6] . Precision measurements of the Z boson decay width [7J exclude models in which the fourth neutrino mass 
eigenstate has a mass less than 45 GeV. If the neutrino masses are very large, however, a fourth generation of fermions 
is not excluded. 

One consequence of the extra fermions is that the ggH coupling is enhanced by a factor of roughly three relative 
to the SM coupling [1, H, 0. The reason for this is that the lowest-order ggH coupling arises from a quark loop. The 
top quark contribution is the largest due to its large coupling with the Higgs boson. In the limit m g 4 3> ran where 
m g 4 is the fourth-generation quark mass, the Higgs boson coupling cancels the mass dependence for each of the three 
propagators in the loop, and the contribution to the ggH coupling becomes asymptotically independent of the masses 
of the two fourth-generation quarks. Each additional fourth-generation quark then contributes as much as the top 
quark, and the ggH coupling is thus enhanced by a factor K e of approximately three. 

The production cross section will be enhanced by a factor of K%. For mj? near the low end of our search range, 
rriH ~ 110 GeV, the gg —> H production cross section is enhanced by roughly a factor of nine relative to the SM 
prediction. This factor drops to approximately 7.5 near the upper end of the search range, mn ~ 300 GeV, assuming 
asymptotically large masses for the fourth-generation quarks. The reason for this drop is that the denominator of 
the enhancement factor, the SM cross section, has a larger contribution from the SM top quark as run nears 2m t . 
The partial decay width for H — »■ gg is enhanced by the same factor as the production cross section. However, 
because the decay H — > gg is loop-mediated, the H — > W + W~ decay continues to dominate for Higgs boson masses 
m H > 135 GeV. 

We consider two scenarios for the masses of the fourth-generation fermions. In the first scenario, the "low-mass" 
scenario, we set the mass of the fourth-generation neutrino m V 4 = 80 GeV, and the mass of the fourth-generation 
charged lepton = 100 GeV, in order to evade experimental constraints [l0| and to have the maximum impact on 
the Higgs boson decay branching ratios. In the second scenario, the "high-mass" scenario, we set m v ± = mp^ = 1 TeV, 
so that the fourth-generation leptons do not affect the decay branching ratios of the Higgs boson. In both scenarios, 
we choose the masses of the quarks to be those of the second scenario in Ref. [9], that is, we set the mass of the 
fourth-generation down-type quark to be 17144 = 400 GeV and the mass of the fourth-generation up-type quark to be 
to„4 = 777^4 + 50 GeV + 101og(m^/115 GeV) GeV. The other mass spectrum of Ref. § chooses = 300 GeV, 
resulting in slightly larger predictions for a(gg —y H) and nearly identical decay branching ratios. We use the next- 
to-next- to-leading order (NNLO) production cross section calculation of Ref. @, which builds on the NNLO SM 
calculations of Refs. p4H18| . 

The CDF and DO Collaborations have searched for the SM Higgs boson in the decays H — )• W + W~ and H — >• ZZ 
using all SM production processes: gg — > H, qq — > WH, qq — > ZH, and vector-boson fusion (VBF) [1914231 ]. The 
results of these searches for the SM Higgs boson cannot be used directly to constrain fourth-generation models, as the 
ggH coupling is enhanced but the WWH and ZZH couplings are not, and the signal acceptances and the backgrounds 
in the multiple analysis channels differ for the various production modes. 

Previously, the CDF and DO Collaborations have published searches for the process gg — > H —y W + W~ , also 
neglecting the WH, ZH, and VBF signal contributions [24-26]. This paper is an update of Ref. [2(|, which provides 
a fourth-generation interpretation of the gg — > H — > W + W~ searches, and a model-independent limit on o_[gg — > 
H — > W + W~). Here we update the searches detailed in Ref. (26[ with those using 8.2 fb^ 1 from CDF [l9f and 
8.1 fb _1 from DO [2l| and include additional channels [2(| [H, [23[ . Interpretations of SM Higgs boson searches within 
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the context of a fourth generation of fermions have been performed by ATLAS [13] and CMS [Hj], which exclude 
140 < m H < 185 GeV and 144 < m H < 207 GeV, respectively. 

We present new limits on a(gg — > H) x B(H — > W + W~) in which the gg —> H production mechanism is considered 
as the unique signal source. These limits are compared to models for Higgs boson production in which the ggH 
coupling is enhanced by the presence of a single additional generation of fermions. In this comparison, the decay 
branching ratios of the Higgs boson are also modified to reflect changes due to the fourth generation relative to the 
SM prediction. While the decays of the heavy quarks and leptons may include W bosons in the final state, we do not 
include these as additional sources of signal. The branching ratios for H — > W + W~ are calculated using HDECAY [29[ 
modified to include fourth-generation fermions The modified Higgs branching ratio to W + W~ is multiplied by 
the cross section [9j to predict the fourth-generation enhanced gg — > H — > W + W~ production rate. We include as 
an approximation to the uncertainties on the decay branching ratios the relative uncertainties predicted for the SM 
branching ratios due to uncertainties in m&, m c , and a s presented in Ref. [30j | . 

The branching ratio for the decay H — > ZZ is also considerable at large masses, having a peak at around m# = 
135 GeV and rising above run — 180 GeV. In order to include a search for gg — > H — > ZZ , we introduce the 
assumption that B(H — > W + W~)/B(H — > ZZ) remains as predicted by the SM, even though the production cross 
section is enhanced. This assumption holds true in fourth-generation models as predicted by the fourth-generation 
version of HDECAY mentioned above. In setting limits on a(gg — > H) x B(H — > W + W~), we include the process 
gg — >• H — >■ ZZ assuming that its signal yield scales in the same way as the W + W~ channel. 

The event selections are similar for the corresponding CDF and DO analyses. Both collaborations select events with 
large $ T and two oppositely charged, isolated leptons, targeting the H — > W + W~ signal in which both W bosons 
decay leptonically. DO selects events containing electrons and/or muons, dividing the data sample into three final 
states: e + e~, e^/x^, and ^ ' fi~ . Each final state is further subdivided according to the number of jets in the event: 
0, 1, or 2 or more ("2+") jets. However, in the fourth-generation interpretation and the model-independent gg — > H 
limits, due to the low S/B in the 2-jet multiplicity bin, only the and 1-jet multiplicity bins were taken into account. 
Decays involving tau leptons are included in two orthogonal ways. A dedicated analysis {[iThad) using 7.3 fb _1 of 
data studying the final state involving a muon and a hadronic tau decay, is included [22J . Final states involving other 
tau decays and mis-identified hadronic tau decays are included in the e + e~, e^/j^, and final state analyses. 

The CDF analysis separates opposite-sign candidate events into five non-overlapping channels. Events are classified 
by their jet multiplicity (0, 1, or > 2), and the and 1 jet channels are further divided according to whether both 
leptons are in the central part of the detector or if either lepton is in the forward part of the detector. The presence of 
neutrinos in the final state prevents event-by-event reconstruction of the Higgs boson mass and thus other variables 
are used for separating the signal from the background. For example, the angle A(f>(£, £) in signal events is smaller 
on average than that in background events, the missing transverse momentum is larger, and the total transverse 
energy of the jets is lower. The CDF analyses use neural-network (NN) outputs as the final discriminants; the list 
of network inputs includes likelihood ratio discriminant variables constructed from matrix-element probabilities [2jJ . 
The DO e + e~ , e ± ^ T , and final state channels use random forest boosted decision tree (BDT) outputs as the 

final discriminants while the fiThad channel uses neural networks. Both experiments use separate NNs and BDTs for 
the different jet multiplicities to distinguish the gg — > H signal from the backgrounds for each of the test masses, 
which are separated by increments of 5 or 10 GeV. 

CDF introduces a search for gg — > H — > ZZ — > £ + £~£' + £'~ using 8.2 fb _1 of data [20| . The discriminant variable is 
the reconstructed mass of the four leptons muu- The version of this search used for the Standard Model interpretation 
also includes the expected yield for W H, ZH, and VBF production, but these production modes are assumed not 
to contribute when interpreting the results in this paper. 

DO also includes for the first time channels in which one of the W bosons in the H — > W + W~ process decays 
leptonically and the other decays hadronically (23[. Electron and muon final states are studied separately, each with 
5.3 fb _1 of data. Random forests are used for the final discriminants. 

Both CDF and DO conduct all searches included in this paper in the range 110 < mn < 300 GeV. The details of the 
signal and background estimations and the systematic uncertainties are provided in Refs. [19- 23]. We set limits on 
a (99 ~ > H ) x B(H —> W + W~) as a function of m#. We use the same two statistical methods employed in Ref. [3lj |. 
namely the modified frequentist (CL S ) and Bayesian techniques, in order to study the consistency of the results. 
Each method is applied at each test mass to calculate an observed upper limit on cr(gg —> H) x B(H — > W + W~). 
Pseudo-experiments drawn from systematically varied background-only predictions are used to compute the limits we 
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expect to obtain in the absence of a signal. 

Correlated systematic uncertainties arc treated in the same way as they are in Ref. [3l|. The sources of correlated 
uncertainty between CDF and DO are the total inelastic pp cross section used in the luminosity measurement, the SM 
diboson background production cross sections (WW, WZ, and ZZ), and the tt and single top quark production cross 
sections. Instrumental effects such as trigger efficiencies, lepton identification efficiencies and misidentification rates, 
and the jet energy scales used by CDF and DO remain uncorrelated. To minimize the degrading effects of systematics 
on the search sensitivity, the signal and background contributions are fit to the data observations by maximizing a 
likelihood function over the systematic uncertainties for both the background-only and signal+background hypothe- 
ses [Hj]. When setting limits on cr(gg — > H) x B(H — > W + W~), we do not include the theoretical uncertainty on 
the prediction of cr(gg — > H) x B(H — > W + W~) in the fourth-generation models since these limits are independent 
of the predictions. When setting limits on ran in the context of fourth-generation models, however, we include the 
uncertainties on the theoretical predictions as described below. 

Since both the CDF and the DO analyses categorize events into separate channels based on the number of recon- 
structed jets, we include the systematic uncertainties on each jet category according to the procedure recommended 
by Ref. [331 ] . This procedure treats the factorization and renormalization scale uncertainties as independent in the 
inclusive gg — »■ H, gg — > H + 1 or more jets, and gg — > H + 2 or more jets calculations. We propagate these un- 
certainties through to the calculations of the exclusive H + jets, H + 1 jet, and H + 2 or more jets calculations 
required to predict the yields in the fourth-generation models. When setting a limit on a(gg — > H — y W + W~), we 
do not include a theoretical uncertainty on the prediction of the quantity that we are constraining experimentally. 
Nonetheless, we do require theoretical input for the differential distributions of the number of jets in gg — > H pro- 
duction. We therefore use only the factorization and renormalization scale uncertainties on the H + I or more jets 
and H + 2 or more jets inclusive calculations when setting a limit on the inclusive process. We follow the PDF4LHC 
prescription [34| for evaluating the parton distribution function (PDF) uncertainties on the signal predictions. QCD 
scale uncertainties that affect the PDF predictions are included as part of the total scale uncertainty, added linearly 
to the other components of the scale uncertainty. The remaining component of the PDF uncertainty is treated as 
uncorrelated with the QCD scale uncertainty. 

The scale choice affects the px spectrum of the Higgs boson when produced in gluon-gluon fusion, and this effect 
changes the acceptance of the selection requirements and also the shapes of the distributions of the final discriminants. 
The effect of the acceptance change is included in the calculations of Ref. [35] and Ref. [36J , as the experimental require- 
ments are simulated in these calculations. The effects on the final discriminant shapes are obtained by reweighting the 
Pt spectrum of the Higgs boson production in our Monte Carlo simulation to higher-order calculations. The Monte 
Carlo signal simulation used by CDF and DO is provided by the LO generator pythia 37f which includes a parton 
shower and fragmentation and hadronization models. We reweight the Higgs boson px spectra in our pythia Monte 
Carlo samples to that predicted by hqt [38| when making predictions of differential distributions of gg — » H signal 
events. To evaluate the impact of the scale uncertainty on our differential spectra, we use the RESBOS [39| generator, 
and apply the scale-dependent differences in the Higgs boson pt spectrum to the hqt prediction, and propagate these 
to our final discriminants as a systematic uncertainty on the shape, which is included in the calculation of the limits. 

The combined limits on a(gg — > H) x B(H — y W + W~) are listed in Table [TT1 for both the CL S and the Bayesian 
methods, and are shown in Fig. [TJ along with the fourth-generation theory predictions for the high-mass and low-mass 
scenarios. The uncertainty bands shown on the low-mass theoretical prediction are the quadrature sum of the MSTW 
2008 40] 90% C.L. PDF uncertainties and the factorization and renormalization scale uncertainties from Table I of 
Ref. [9j. The limits calculated with the two statistical methods agree within 10% for all Higgs boson mass hypotheses. 
An excess is seen in the neighborhood of m# = 140 GeV. This excess is localized to the DO H — y W + W~ — > e + ^ e /i~P M 
search (including charge conjugate final states) [2l[. The excess is not as pronounced in the SM-optimized search 
from DO [2l[. Combination with CDF, which sees no significant excess at ttih = 140 GeV results in a total local 
significance somewhat in excess of two standard deviations. The trials factor, or look-elsewhere effect, further dilutes 
the significance, in this case by a large factor due to the large range of ran explored. We do not claim that this excess 
is statistically significant. 

In order to set limits on ran in these two scenarios, we perform a second combination, including the uncertainties on 
the theoretical predictions of cr(gg —> H) x B(H — >■ W + W~) due to scale and PDF uncertainties at each value of ran 
tested. The resulting limits are computed relative to the model prediction, and are shown in Fig. [5] for the low-mass 
scenario, which gives the smaller excluded range of m#. In this scenario, we exclude at the 95% C.L. a SM-like 
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FIG. 1: The Tevatron combined observed (solid black lines) and median expected (dashed black lines) 95% C.L. upper limits 
on o(gg — > H) x B(H — > W + W~). The shaded bands indicate the ±1 standard deviation (s.d.) and ±2 s.d. intervals on 
the distribution of the limits that are expected if a Higgs boson signal is not present. Also shown the prediction for a fourth- 
generation model in the low-mass and high-mass scenarios, 4G (Low mass) and 4G (High mass) respectively. The hatched 
areas indicate the theoretical uncertainty from PDF and scale uncertainties. The lighter curves show the high-mass theoretical 
prediction. 



Higgs boson with a mass in the range 124 - 286 GeV. Using the median limits on <r(gg — > H) x B(H — > VF+iy - ), 
expected in the absence of a signal to quantify the sensitivity, we expect to exclude the mass range 120 - 267 GeV. 
In the high-mass scenario, which predicts a larger B(H — )• W + W~) at high mjj than that predicted in the low-mass 
scenario, we exclude at the 95% C.L. the mass range 124 - 300 GeV and expect to exclude the mass range 120 - 
290 GeV. The upper edge of the observed exclusion range in the high-mass scenario is determined by the fact that 
we did not perform the search for Higgs bosons with mass exceeding 300 GeV. 

In summary, we present a combination of CDF and DO searches for the gg H — > W + W~ and gg — > H — > ZZ 
processes and set an upper limit on a(gg — > H) x B(H —> W + W~) as a function of rn#, assuming the standard 
model ratio of B(H — > W + W~)/B(H — > ZZ). We compare these limits with the prediction of the minimal SM with 
a sequential fourth generation of heavy fermions added on, and exclude at the 95% C.L. the Higgs boson mass range 
124 < m H < 286 GeV, with an expected exclusion of 120 - 267 GeV. 
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FIG. 2: The Tevatron combined observed (solid black lines) and median expected (dashed black lines) 95% C.L. upper limits 
on the ratio of the Higgs boson production cross section to that predicted by the 4G (Low mass) model; the uncertainties in 
the signal prediction are included in the limit. Also shown is the prediction of the signal rate in the high-mass scenario, divided 
by that of the low-mass scenario. 
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TABLE I: The fourth-generation enhanced a(gg — > H) listed in fb for the low-mass scenario described in the text with 
m d4 = 300 GeV and m d4 = 400 GeV along with the B(H ->■ W + W~) for Higgs masses between 110 - 300 GeV. The 
uncertainty on the predicted cross section from variations in the PDF and factorization and renormalization scale are also 
listed in percentage where these have been determined from the MSTW 2008 90% C.L. uncertainty and by modification of the 
scale by factors of 1/2 and 2, respectively. The inclusive uncertainties are listed here, but the differential uncertainties in the 
jet categories are used in the calculation of the limits. 

m H cr(gg — >• H) a(gg -> H) uncert. uncert. uncert. uncert. BR(H -> W + W~) BR(H -»■ W + W~) 

(GeV) m d4 = 300 GeV m d4 = 400 GeV PDF up(%) PDF down(%) (j, up(%) (j, down(%) m d4 = 300 GeV m d4 = 400 GeV 
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5094.6 


5050.4 


13 


-12 


12 -8 


0.44 


0.44 


150 


4540.5 


4498.5 


13 


-12 


12 -8 


0.55 


0.55 


155 


4055.6 


4017.6 


13 


-12 


12 -8 


0.68 


0.68 


160 


3630.2 


3595.1 


13 


-13 


12 -8 


0.85 


0.85 


165 


3253.7 


3220.7 


14 


-13 


12 -8 


0.91 


0.91 


170 


2924.1 


2893.2 


14 


-13 


12 -8 


0.89 


0.88 


175 


2633.9 


2604.4 


14 


-13 


12 -8 


0.86 


0.86 


180 


2376.7 


2348.9 


14 


-13 


12 -8 


0.83 


0.83 


185 


2147.2 


2121.5 


15 


-13 


12 -8 


0.74 


0.74 


190 


1943.9 


1919.7 


15 


-14 


12 -8 


0.69 


0.69 


195 


1763.2 


1740.2 


15 


-14 


12 -8 


0.66 


0.66 


200 


1601.8 


1580.0 


15 


-14 


12 -8 


0.65 


0.65 


210 


1328.1 


1308.4 


16 


-14 


12 -8 


0.62 


0.62 


220 


1107.7 


1089.6 


16 


-15 


12 -8 


0.60 


0.60 


230 


928.61 


912.21 


17 


-15 


12 -8 


0.59 


0.59 


240 


782.52 


767.44 


17 


-15 


12 -8 


0.58 


0.58 


250 


662.60 


648.81 


18 


-16 


12 -8 


0.58 


0.58 


260 


563.53 


550.90 


19 


-16 


12 -8 


0.57 


0.58 


270 


481.49 


469.93 


19 


-16 


12 -8 


0.57 


0.57 


280 


413.24 


402.68 


20 


-17 


12 -8 


0.58 


0.58 


290 


356.39 


346.53 


21 


-17 


12 -8 


0.58 


0.58 


300 


308.70 


299.71 


21 


-17 


12 -8 


0.58 


0.58 
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TABLE II: The observed and median expected 95% C.L. upper limits on a(gg — >• H)xB(H — > W + W~) for inn between 110 GeV 
and 300 GeV, obtained with the Bayesian and CL S methods. Also listed are the a(gg — > H) x B(H — > W + W~) predictions 
of the low-mass and the high-mass fourth-generation scenarios discussed in the text for a fourth-generation down-type quark 
mass of 400 GeV. All limits and predictions are presented in pb. 
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